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The superfluid ®He formed by spin-triplet p-wave Cooper pairs is a typical topological superfluid. 

In the snperfluid ®He B-phase, several kinds of vortices classified by spatial symmetries Pi, P 2 , and 
P 3 are produced, where Pi is inversion symmetry, P 2 is magnetic reflection symmetry, and P 3 is 
magnetic yr-rotation symmetry. We have calculated the vortex bound states by the Bogoliubov-de 
Gennes theory and the quasiclassical Eilenberger theory, and also clarified symmetry protection of 
the low energy excitations by the spatial symmetries. On the symmetry protection, P 3 symmetry 
plays a key role which gives two-fold degenerate Majorana zero modes. Then, the bound states in 
the most symmetric o vortex with Pi, P 2 , and P 3 symmetries and in w vortex with P 3 symmetry 
have the symmetry protected degenerate Majorana zero modes. On the other hand, zero energy 
modes in v vortex, which is believed to be realized in the actual B-phase, are not protected, and 
in consequence become gapped by breaking axial symmetry. The excitation gap may have been 
observed as the variation of critical velocity. We have also suggested an experimental setup to 
create o vortex with Majorana zero modes by a confinement and a magnetic field. 

PACS numbers: 67.30.he 


I. INTRODUCTION 

Superfluid ®He is a condensate of spin-triplet p- 
wave Cooper pairs, which undergoes complex symmetry 
breaking, in addition to U(\)n gauge symmetry break¬ 
ing Q. Without magnetic fields, the superfluid phases 
of ®He consist of the nodal gapped A-phase in a high- 
temperature and high-pressure region and the B-phase 
with an isotropic gap in another wide parameter region. 
In the full gapped B-phase, the residual symmetry group 
is 

= SO{J.)l+s xTxC, (1) 

where iS'0(3)i,+s: is simultaneous rotation symmetry of 
orbital and spin spaces, T is time-reversal symmetry, and 
C is particle-hole symmetry. Three-dimensional spin- 
triplet superfiuids with T and C, such as the B-phase, are 
topological phases belonging to class Dill in the Altland- 
Zrinbauer symmetry classes In consequence of the 
topological order, the Majorana bound states and heli¬ 
cal spin current exist on surfaces of the superfluid ®He 
B-phase [l-@. 

The topological phases are classified by whether there 
are discrete symmetries T and C 0. It has been known 
that also spatial discrete symmetries in point group can 
give rise to a nontrivial topology of excitations in the 
topological phases 0 - The topological excitations with¬ 
out a gap are robust if the symmetries are preserved. 
Host matters of the symmetry protected excitations are 
called topological crystalline insulators 0 or topological 
crystalline superconductors The symmetry pro¬ 

tected excitations are mainly discussed on gap less bound¬ 
ary states with reflection symmetries [3, . Recently, 


the topological classification is generalized for excitations 
in a defect with a two-fold symmetry including magnetic 
point group symmetry 

An example of the symmetry protected excitation 
in the superfluid ®He B-phase is the surface Majorana 
bound state under magnetic fields [l^. The gapless Ma¬ 
jorana bound state is preserved under weak magnetic 
fields parallel to the surface although time-reversal sym¬ 
metry is broken [l 3 , [3 • This system has magnetic tt- 
rotation symmetry instead of time-reversal symmetry be¬ 
cause the flipped magnetic field by time-reversal opera¬ 
tion is recovered by the 7r-rotation around the axis per¬ 
pendicular to the surface. The gapless Majorana bound 
state under magnetic fields is due to the protection by 
the magnetic 7r-rotation symmetry. 

The symmetry protected Majorana excitations are also 
expected in the vortex bound state, where time-reversal 
symmetry is broken by the phase winding of the vortex. 
For an axisymmetric vortex in the B-phase, the residual 
symmetry group is 

Hvortex — C/(1 )q xPxC, (2) 

where the generator of the U{1)q symmetry is Q = 
Lz — N/ 2 for a singly quantized vortex along the z- 
axis [3, . The residual spatial symmetry P is P 2 X P 3 

for the most symmetric vortex called o vortex , where 
P 2 is magnetic reflection symmetry on a plane includ¬ 
ing the vortex line and P 3 is magnetic 7r-rotation sym¬ 
metry around a axis perpendicular to the vortex line. 
The combination of P 2 and P 3 gives inversion symmetry 
Pi = P 2 P‘i- A vortex with a symmetry P = Pi, P 2 , or 
P 3 is called u, v, or w vortex, respectively, and the lowest 
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symmetric vortex without the symmetry P is called uvw 
vortex. 

In experiment, two types of vortices, VI vortex in a 
high-temperature and high-pressure region and V 2 vor¬ 
tex in a low-temperature and low-pressure region, were 
observed within the B-phase [2^. The transition line 
between the vortex states was determined by the NMR 
measurement [ 22 l - [^ and the measurement of critical ve¬ 
locity [2^ [2^ . Since the transition pressure on the super¬ 
fluid critical temperature is reproduced by the Ginzburg- 
Landau theory, VI vortex and V 2 vortex are believed to 
be axisymmetric v vortex and non-axisymmetric v vor¬ 
tex, respectively [ 13 , Breaking of the axial symmetry 
over the transition from VI vortex to V 2 vortex was con¬ 
firmed by the observation of the Goldstone mode related 
to twisting of the anisotropic vortex core [ 1 ^. The exis¬ 
tence of P2 symmetry which is a direct evidence to form 

V vortex, however, has not been verihed. 

For the identification of the kind of vortex in the super¬ 
fluid ^He B-phase, the symmetry protected vortex bound 
state can be utilized. Indeed, we will demonstrate that P3 
symmetry guarantees the existence of Majorana zero en¬ 
ergy modes in the vortex bound state. On the other hand, 
P2 symmetry does not provide symmetry protected exci¬ 
tations in a vortex. Therefore, many zero energy modes 
in axisymmetric v vortex [s^ are gapped out by the ax¬ 
ial symmetry breaking. The formation of the gap in low 
energy excitations through the vortex transition implies 
the realization of v vortex with only P2 symmetry. 

This paper is arranged as follows: In Sec. II, we for¬ 
mulate the Bogoliubov-de Gennes (BdG) theory and the 
quasiclassical theory. Vortex bound states are calculated 
by the BdG theory with order parameters (OPs) which 
are self-consistently obtained by the quasiclassical the¬ 
ory. The possible spatial symmetries Pi, P2, and P3 for 
an axisymmetric vortex in the B-phase are summarized 
in Sec. III. The spatial symmetries play an important role 
in the topological classification of the vortex bound state. 
In Sec. IV, we discuss the bound state in the most sym¬ 
metric o vortex with Pi, P2, and P3 symmetries. The vor¬ 
tex bound state has symmetry protected Majorana zero 
modes; however, o vortex is not realized in the actual B- 
phase. The v vortex with P2 symmetry is believed to be 
realized in the B-phase, whose bound state is discussed 
in Sec. V. There is a difference between low energy exci¬ 
tations in axisymmetric v vortex and non-axisymmetric 

V vortex in consequence of that P2 symmetry does not 
protect zero energy modes. We also discuss the bound 
state in w vortex with P3 symmetry in Sec. VI. Gon- 
cerning the OP, the difference between w vortex and v 
vortex is only a phase of induced components of the OP 
around the vortex. However, the bound state in w vor¬ 
tex dramatically changes from that in v vortex and has 
symmetry protected Majorana zero modes owing to P3 
symmetry. We devote the final section to the summary in 
which we also mention topologically trivial bound states 
in uvw vortex and u vortex. 


II. FORMULATION 
A. Bogoliubov-de Gennes theory 

We have numerically calculated the vortex bound 
states for some kinds of vortices in the B-phase by the 
BdG theory. The BdG theory gives discretized modes 
in the vortex bound state at A^/Pp intervals, where 
A is a superfluid gap and Pp is the Fermi energy [ 3 l|. 
Note that the discreteness can be negligible when we con¬ 
sider physical quantities in the superfluid ^He owing to 
A/Pp ^ 10 “^ [I|; however, we use the BdG theory in or¬ 
der to discuss the symmetry protection for the discretized 
modes in the vortex bound state. 

The BdG equation in spin-triplet superfluid states is 
given as [ 13 ] 



f Hri,r2) 
V-A't'(ri,r2) 


Hri,r2) \ 
-/i'^(ri,r2)y 


Uu{r2) 


- (T *3 ) , 


( 3 ) 


where the quasiparticle wave function with the v- 
th excited state with the eigenvalue E^, is u^ir) = 
[ul(r),ui{r), uj(r),r!;j;(r)]"^. The single particle Hamil¬ 
tonian and the OP are described as 


Hri,r2) 


A(ri,r2) 


2 m 


E-p 


Siri - r2)i. 


A^|(ri,r2)\ 

VAit(’’i>’’2) A;|(ri,r2)y’ 


( 4 ) 

( 5 ) 


where m is mass of the condensed particle and 
Aaa'iri,r2) = V{r')Y,^u^{ri)[v^' {r2)]* f{E^) with the 
interparticle interaction V(r') at r' = \ri — r2\ and the 
Fermi distribution function /(Pj^). 

Here, we consider the bound state in a vortex along the 
2-axis; then the wave number becomes a well-defined 
quantum number. The BdG equation for the A:^-resolved 
two-dimensional ( 2 D) form is given as [sH ] 



(PiJ P2) 

-A^. (Pi,P2) 


Afej {Pi, P2) 


Uu,kAP2) 


— (Pi)j (6) 


where 


^kAPi^ P2) — 


2 m ^ ^ 


S{p,-pAl ( 7 ) 


with V2P) = + 9 y_^. The 2 D form of the Fermi en¬ 

ergy E^^^kz) = (/i^/2m)(fc| — kl) reflects the fc^-cross 
section of the Fermi surface, where fcp is the Fermi wave 
number. Although the OP Ak^iPi, P2) should be calcu¬ 
lated self-consistently with the quasiparticle wave func¬ 
tion (pi) and the eigenvalue Pj/,fc^, an approximate 
solution can be derived from the quasiclassical theory for 
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superfluid states with a small gap A/Ep ^ 1. We ex¬ 
pand the OP to the Fourier integral with the relative 
coordinate p' = pi — p 2 as 

- f 2n / - 

A,APi,P 2 ) = J ■^A(fe,p)r(fc), ( 8 ) 

where p = {pi + P 2)/2 is the center-of-mass coordinate 
and = {kx,ky) is the 2D component of the relative 
momentum k. In this model, the Gaussian factor T{k) = 
g-(,K -KpKp indicates that the pairing interaction is non¬ 
zero near the Fermi surface k = kp with the range of the 
interaction [s^. The OP A{k,p) is obtained by the 
quasiclassical theory in Sec. II B. 

If the fc 2 -resolved quasiparticle wave function Uu,k^ = 

satisfies the condition 

= {'^1,-k^) and , the creation 

and annihilation operators are equivalent for the zero 
energy quasiparticle with k^ = 0 [3^ . Then, the 

quasiparticle is a Majorana zero mode. 


B. Quasiclassical theory 

The spatial structure of OP with a vortex has been 
calculated by the quasiclassical theory, which is valid for 
superfluids and superconductors with A ^ Ep, such as 
the superfluid ^He. We have found self-consistent solu¬ 
tions of the OP A{k,p) with the quasiclassical Green’s 
function g{k, p,ujn) by the Eilenberger equation [Sbl - I^ 


is more convenient than the description by the d-vector 
for the OP with an axisymmetric vortex, namely, 


A(fc,p) 


(-V2C+{k,p) Coikp) \ 

{ Coik,p) V2C.{k,p)J^ 


( 12 ) 


where C± = {dx T 'idy)/'/2 and Co = dz- Each coeffi¬ 
cient can be expanded in projections of orbital angular 
momentum. 


C'a(fe, p) — Ca_|_(p)fc_|_ -|- Coo(p)^0 + C'a-(p)fc-, (13) 

with a = 0 or ±, where k± = (kx ± iky)/\/2 and fco = 

kz El. 

The self-consistent condition for A is given as 

A{kp)=No7rkBT ^ (V{k,k')f{k',p,ujr,))_,, 

(14) 

where Nq is the density of states in the normal state, uJc is 
a cutoff energy setting uic = ^OksTc with the critical tem¬ 
perature Tc, and ■ ■)k indicates the Fermi surface aver¬ 
age. The pairing interaction V{k,k) = 3gk ■ k , where 
g is a coupling constant with the relation (pAq)“^ = 
ln(r/rc) + We solve Eq. © and 

Eq. (fTTl) alternately at T = 0.2Tc, and obtain a self- 
consistent solution. Then, we use the self-consistent OP 
A{k, p) after the replacement of k by k/kp as A{k, p) in 
Eq. 


- ihv{k) ■ Vg{k,p,ujn) 


f fWnl 
\hA{k,p) 


-A(fc,p)\ 
-iUJni ) 


,g{k,p, ujn) 


(9) 


where the wide hat indicates the 4x4 matrix in particle- 
hole and spin spaces. The quasiclassical Green’s function 
is described in particle-hole space by 


g{k,p,u!n) = -ITT 


f 9{k,J>,UJn) 
\-il{k,p,UJn) 


-9{k,P,UJn)J 


where = (2n -I- IjTr/csP is the Matsubara frequency 
and k is the normalized relative momentum on the Fermi 
surface. The quasiclassical Green’s function satisfies the 
normalization condition 'g^ = —tt^I. The Fermi velocity 
is given as v{k) = vpk on a three-dimensional Fermi 
sphere. 

The spin-triplet OP is defined by 

A{k,p) = id{k,p) ■ ady, (11) 

where & is the Pauli matrix. The d-vector is perpendicu¬ 
lar to the spin S' of a Cooper pair, namely, d- S = 0. The 
description using projections of spin angular momentum 


III. SPATIAL SYMMETRIES FOR AN 
AXISYMMETRIC VORTEX 

For the OP with an axisymmetric vortex, the coeffi¬ 
cients in Eq. (US are described by Cab{p) = Cab{p)d'^'^'‘'^, 
where p is a radial distance from a vortex core and (j) is 
an azimuthal angle in the a;?/-plane. Since the OP has 
U{1)q symmetry with Q = Lz + Sz — N/2 as Eq. ([2|), the 
phase winding number nab in each coefficient satisfies 

na-t-l,fe — ‘k^a^b Ij 
^a,fe+l — ‘^a,b 1- 

The following coefficients are finite in the bulk B-phase: 

C+_=Coo = G_+ = Ab, ( 16 ) 

where Ab is the amplitude of a gap in the bulk, because 
the total angular momentum J = i-|-S = Oin the 
B-phase El- With a singly quantized vortex, the OP 
has 27r-phase winding around the vortex, namely, phase 
winding numbers 

n+_ = noo = n—I-= 1 . ( 17 ) 
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Then, the general description of the coefficients with a 
singly quantized axisymmetric vortex in the superfluid 
^He B-phase is 

(C++{p)e-^* C+o(p) C'+_(p)e*n 

C{p) = Co+{p) Cooip)e^^ Co-{p)e^^^ , 

VC_+(p)e*^ C-oipy^* C—{p)e^^y 

(18) 


gives 

P3Cik,p,^) (P*)-' 

-c%{p) Cl_ip)P^ \ /k+\ 

= -cs+ip) cSoip)P^ ho ■ 

VC'l+(p)e*‘^ -C*_o{p)e^^^ C*__{p)e^^*) \k^) 

(23) 


where Cab approaches Ab for a + 5 = 0 and vanishes for 
a + 6 0 when p —?► oo. 

The B-phase with a vortex can have the additional 
symmetry P as Eq. ©. Under the symmetry operation, 
the semi-classical BdG Hamiltonian 


P) 


(/Hk) A{k,p)\ 
\A^k,p) -fP'i-k))’ 


(19) 


satisfies 


IV. o VORTEX 

The o vortex is the most symmetric vortex in the super¬ 
fluid ^He B-phase. Since it has the all possible discrete 
symmetries Pi, P 2 , and P 3 , several coefficients should 
vanish as 

/C++{p)e-^^ 0 C+.{p)P<^\ 

C{p) = 0 Coo(p)e'’' 0 , (24) 

Vc_+(p)e*^ 0 C__(p)eO*V 


VHBdG{k,p)V-^ = HBdG{k',p'), ( 20 ) 

where the symmetry operator is P = diag('P,'P*) and k 
and p are transformed into k' and p', respectively, by 
the operation. The normal state Hamiltonian h{k) = 
{h?/2m){k^ — kp)i implies Ph{k)P~^ = h{k). Since 

the operation PA{k,p)(p*'^ can be regarded as 

PC{k,p){P*) the inversion symmetry operator Pi 
acts on Ca{k, p, cji) as 

PiC{k, p, 4>) {Pir^ = c(-fc, p, </)+^) 

/C++(p)e-*^ -G+o(p) C+_(p)e*<^ \ (k+\ 

= -Co+(p) Coo(p)e*‘^ -Co_(p)e2*<^ h^ ' 

VC_+(p)e*^ -C_o(p)e2*^ ) \k-) 

( 21 ) 

The symmetry P 2 is the combined symmetry with 
the time-reversal and mirror reflection on a plane in¬ 
cluding the vortex line, namely the magnetic reflec¬ 
tion symmetry. The time-reversal operator T acts as 
TCa,b{p,(j)){T*)~^ = -Cl^ _i^{p,(j)) and the mirror re¬ 
flection operator M acts as Vf C'o,&(p, </>) (Af *)~^ = 
-C-a-bip,-4>)l therefore, 

P2C{k,p,4,) {p*)-^ 

(cx+{p)e-^^ c;o(p) c;_(p)e*n /fc+\ 

= CoV(p) C*o(p)e*^ Q_(p)e 2 *^ • 

VC'*+(p)e*'^ C'*o(p)e 2 *‘^ C'*_(p)e3*-^/ \k-) 

( 22 ) 

The symmetry P 3 is the magnetic 7 r-rotation symme¬ 
try around an axis perpendicular to the vortex line. 
The TT-rotation operator P acts as PCa,b{p,4>) = 

C-a,-b{PiT^ — 0)- Then, the combination of T and P 


where all remaining coefficients Cab{p) are real and they 

approach (7+- = Coo = C _= Ab and (7++ = C _= 0 

when p — 00 . The component (7++ (C--) is induced 
around the vortex by the spatial variation of the bulk 
component _(G-y.) [ 2 l|. 

The self-consistently obtained OP by the quasiclassi- 
cal theory is shown in Fig. [IJa). The bulk components 

( 7-1 _, Coo, and C _ y rise up with p-linear from a vortex 

core, where Coo recovers the bulk gap Ab with slightly 

shorter length than (7-|_and C _Rises of C++ and 

C _are p-linear and p-cubic, respectively, which reflect 

phase winding numbers n++ = — 1 and n _= 3. This 

difference makes a little variance between (7-|_and C _|- 

via a coupling between the same orbital state. 

The discretized eigenvalues of the vortex bound state, 
which is derived from the BdG theory with the self- 
consistent OP, are shown in Figs. Hd) and HKe). In 
Fig. HKd), the eigenvalues for k^ = 0 are classified into 
quantized orbital angular momentum 1. The low energy 
eigenvalues are discretized at the order of Ag/Fp inter¬ 
vals, where we take Ab/Ff = 0.1. The vortex bound 
state has spin degenerate exact two zero energy modes 
at Z = 0 indicated by an arrow. The wave functions 
of the degenerate zero energy modes for f-spin and 4 ,- 
spin states are shown in Figs. Hb) and [ijc), respec¬ 
tively. Since the orbital chirality of the t"Spin state 
is antiparallel to the vorticity but the orbital chirality 
of the 4 .-spin state is parallel to the vorticity, uj(p) = 

vlip) oc Jo(fcF/o)exp(-p/^) and u^{p) = v^{p) 


Ji(fcFp) exp(—p/^)e®‘^ [39|, l40|, where Jn is the n-th or¬ 
der Bessel function and {k+^)~^ = Ap/Fp = 0.1. These 
zero energy modes are Majorana zero modes because the 
quasiparticle wave functions for both spin states satisfy 

the condition Mq = I 1 and wg = I Uq 1 , as demon¬ 
strated in Fig. [T](b) and (He). 

The Zc^-dispersion of the vortex bound state for var¬ 
ious Vs is shown in Fig. [He), where two different spin 
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FIG. 1. (Color online) Numerical results for o vortex, (a) OP 

configuration Cab{p) for finite components C-i_, Coo, C_ 

C++, and C_Real part of wave functions ■Uo(p) Vo(p) 

with = 0 for the t-spin state (b) and 4.-spin state (c). 

(d) Eigenvalues for quasiparticles with kz = 0. The degener¬ 
ate eigenstates indicated by an arrow correspond to the two 
wave functions in (b) and (c). (e) fc^-dispersion of the vor¬ 
tex bound state whose color denotes angular momentum I of 
quasiparticles. 


states exist for each 1. As demonstrated in Appendix [Al 
the I = 0 branches linearly cross the zero energy at 
fcj = 0 as E± oc ±kz. Their wave functions are dis¬ 


tinguished into u± = (^Uq ±Uq'^ /v^ + 0 {kz/kF), where 


Uq = 




o') and Un = ^ 0 ,itn ,0 




are 


1 


the wave functions of the Majorana zero modes in the 
f-spin and 4,-spin states, respectively. Note that the in¬ 
version symmetry Pi maps each eigenstate to another 
one, namely, u± = ■piu+. 

The spin degenerate Majorana zero modes are pro¬ 
tected by spatial symmetry. The vortex bound state with 
P 2 and P 3 symmetries is classified by a topological in¬ 
variant Z which is demonstrated in Appendix |B] by using 
Clifford algebras. Note that Pi symmetry does not pro¬ 
vide symmetry protected defect zero modes in general. 

The appropriate topological invariant for o vortex 
is calculated by the one-dimensional winding number, 
which is evaluated from the BdG Hamiltonian at an in¬ 
finite point from a vortex 


HBdoik, 4>) 


( h{k) A(fe,<^) \ 

\A^k,cj)) -h'^{-k))' 


(25) 


where h{k) = {h?/ 2 m){k^ — k^)! and 


A(fe, 0 ) 


Ab / kx iky kz ^ 

kp y kz kx ~\~ ikyJ 


From combination of P 3 and particle-hole sym- 
inetry, the BdG Hamiltoniajn^ for o vortex obeys 
{CV3)HBdG{kx,—ky,—kz-,—(t>){CV3)~^ = —HBdG{k,(l)) 

where V 3 = iaz^lC and C = ry/C are operators for the 
magnetic 7r-rotation and particle-hole symmetry with the 
complex conjugation operator fC. The combined symme¬ 
try defines the chiral symmetry, THBdG{kx^ ky = 0 ,kz = 
0,4i = 0,7r)f“^ = -HBdG{kx,ky = 0,kz = 0,4> = 0,7r) 
with r = TyCTj, in the symmetric space ky = kz = 0 and 
^ = 0 or TT. The chiral symmetry enables us to introduce 
the one-dimensional winding number as [4l|, 1 ^ 


■p—0,7T _ 


47ri 


dkxti' 


TPbug^^^’ fcy — 0, ^2 — 0, *4^ — 0, ky — O^kz — 0, (/> — 0, tt) 


(26) 


For o vortex, the one-dimensional winding number is 
evaluated as = 2 and w'^ = —2. The difference of 
the winding {w^ — w'^) 12 = 2 provides the Z topological 
invariant for o vortex. This topological number corre¬ 
sponds to the number of the zero energy states at kz = 0 . 
(More precisely, the Z topological invariant is equal to 
the index trF of the quasiparticle states at kz = 0. Since 
the chiral symmetry requires that the zero energy states 
at fcz = 0 are eigenstates of F, the index trF reduces 
to the difference between the number of the zero energy 
states with eigenvalue F = -|-1 and that with F = — 1 
Hence, if (w° — w^)/2 = N, there exist at least \N\ zero 
energy states at kz = 0.) The obtained topological num¬ 
ber, i.e. {w° — w'^)/2 = 2 , guarantees the existence of two 
zero energy states at = 0 , which are indeed realized 


as the two I = 0 zero modes in Fig. [TJ 

Note that the second Chern number, which character¬ 
izes bound states in a line defect for the symmetry class 
D, vanishes in the presence of Pi symmetry. The second 
Chern number Ch 2 is obtained by [H, Q 

Ch 2 = ^ y dkd<l)P^hr[UJ,k], (27) 

where + z[a„,a; 3 ]™” is the cur¬ 

vature of the non-Abelian Berry connection Oq with 
a = (^i,i) = {(j},kx,ky,kz), and the non-Abelian Berry 
connection is given by 

=-i{m,k,(j)\da\n,k,(l)), (28) 
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with eigenstates |to, fe, </>) and \n, k, (p) of the BdG Hamil¬ 
tonian HBdcik, 4>)- The indices m and n label quasi¬ 
particle states with negative energies. For o vortex, the 
eigenstates have Pi symmetry which gives the periodic¬ 
ity of Berry curvature, f^i{k, (/> + tt) = —f^i{—k, p) and 
fjk{k,(j) + TT) = fjk{—k,cj)). Therefore, the integral in 
Eq. (l27l) yields Ch 2 = 0. Since the second Chern num¬ 
ber generally vanishes when there is Pi symmetry, it is 
not appropriate for the topological number of the bound 
states in o vortex. 


V. V VORTEX 

In the actual B-phase, it is believed that two kinds of 
V vortices are realized. One v vortex has axial symmetry 
and the other v vortex breaks the axial symmetry. In 
this section, we show the difference of the vortex bound 
states between the two kinds of v vortices. 


A. Axisymmetric v vortex 

The axisymmetric v vortex has the magnetic reflection 
symmetry P 2 . From Eq. (1^ . all coefficients Cab{p) in 

Eq. (fTSll are real, where = Coo = C_= Ab and 

the other coefficients vanish when p —> 00 . Since Co+ and 
C+o corresponding to the A-phase and ,5-phase compo¬ 
nents, respectively, do not have any phase winding, they 
can compensate a vortex core. 

The self-consistently obtained OP by the quasiclassi- 
cal theory is shown in Fig. ma). The A-phase component 
Co+ compensates the vortex core with larger amplitude 
than that of the bulk components in the B-phase. The ,5- 
phase component C+o also compensates the vortex core; 
however, the amplitude is smaller than that of Co+ and 
the sign is opposite. Away from the vortex core, Co+ 
and C+o are identical to Co- and C_o, respectively, in 
order that total angular momentum becomes zero. A 
core radius of the v vortex, which is characterized by 
the healing length of the bulk components, is longer than 
that of o vortex because a loss of the condensation en¬ 
ergy at the V vortex core is small. Other and C_ 

components are also induced slightly, which change the 
sign away from the vortex core. 

The discretized eigenvalues derived from the BdG 
equation by using the self-consistent OP are shown in 
Figs, ^d) and H^e). In Fig. [DJd), the eigenvalues for 
kz = 0 are classified into quantized orbital angular mo¬ 
mentum I along the vortex line. The quasiparticles with 
I = 0 and kz = 0 have the finite energy by the 7r/2-phase 
shift owing to the A-phase component compensating the 
vortex core [see Appendix!^. Discretized eigenvalues for 
kz = 0 rnay be situated on the zero energy with finite 
±/, accidentally, according to a core radius. The wave 
function of an accidental zero energy mode, indicated by 
an arrow in Fig. Hd), is shown in Figs. [H[b) and (He) 



FIG. 2. (Color online) Numerical results for axisymmetric v 
vortex, (a) OP configuration Cab{p) whose legends are identi¬ 
cal to that in Fig. UK a), (b) and (c) Real part of wave function 
uZ,k^=o{p) a-nd vZ^k^=oip) for a state with « 0. (d) 

Eigenvalues for quasiparticles with kz = 0. The eigenstate in¬ 
dicated by an arrow corresponds to the wave function in (b) 
and (c). (e) fc 2 -dispersion of quasiparticles with I > 0 whose 
color denotes angular momentum I of quasiparticles. 


for f-spin and ),-spin parts, respectively. The quasipar¬ 
ticle wave function satisfies ul ^^=0 ~ kz=o) while 

k =0 ^ (^'^tkz=o) ’ therefore, the quasiparticle is not 

the Majorana zero mode [H, [s^. It is a natural conse¬ 
quence because the particle-hole operation changes the 
sign of orbital angular momentum. 

The /cj;-dispersion of the vortex bound state for ^ > 0 
is shown in Fig. He). Eigenvalues for I' = —I < 0 are 
given by E' = —E with = —kz owing to the particle- 
hole symmetry. The eigenvalues with I = 0 approach 
the zero energy at kz = Ekp owing to an induced C_o 
component around the vortex, which is different from the 
result by Silaev [see Appendix [C]. Since momentum 
kz is a continuous quantity, the eigenvalues with |/| ;< 
^ckp cross the zero energy at finite kz- The zero energy 
modes with finite kz are, however, also not Majorana 
zero modes because signs of I and kz for the zero energy 
modes are changed by the particle-hole operation, that 
is, the zero energy modes are mapped to different zero 
energy modes by the particle-hole operation. 

The trivial vortex bound state without Majorana zero 
modes can be understood by the topological arguments 
in Appendix |B] The vortex bound states without any 
symmetry protections are clearly shown as gaps of low 
energy excitations in non-axisymmetric v vortex which 
breaks axial symmetry but keeps P 2 symmetry. 



















7 



FIG. 3. (Color online) Numerical results for non- 
axisymmetric v vortex. 2D profile for each OP compo¬ 
nent Cab(p) (a) and the root-mean-square value of the gap, 

|A| = (Tr[AAt])fc/2, (b). (c) Eigenvalues for quasiparti¬ 
cles with kz = 0. (d) fcj-dispersion of the vortex bound state 
whose color denotes angular momentum I of quasiparticles. 

B. Non-axisymmetric v vortex 

The self-consistently obtained OP for non- 
axisymmetric V vortex is shown in Figs. [31(a) and 
ISKb). Far away from the vortex core, non-axisymmetric 
V vortex is aiso described by the coefficients in Eq. (HI; 
however, the originai 47r-phase singuiarity of C-o and 
Co- is spiit into two singuiarities with the 27r-phase 
winding as shown in Fig. |2Ka). This deformation 
breaks 17(1) q symmetry but keeps P 2 symmetry. In 
Fig. |2Ib), we show the root-mean-square vaiue of the 

gap, |A| = (Tr[AA'i])^/2. Since four components 

Co+, C+o, C-o, and Co- are finite at the vortex center, 
substantiai gap opens on the vortex center. In return 
for the gap on the vortex center, finite minima of |A| 
appear on the phase singuiarity of C-o and Co-. 

The discretized eigenvaiues derived from the BdG 
equation by using the seif-consistent OP are shown in 
Figs. Elc) and ISKd). In Fig. ISKc), the eigenvaiues for 
kz = 0 are ciassified into orbitai anguiar momentum I 
aiong the vortex iine. The /c^-dispersion of the vortex 
bound state with vaiues of I is shown in Fig. EKd). Note 
that I is not a quantum number but it is caicuiated by 

l = -ihj dpulj^^ip)d^u^,kAp), (29) 

for each eigenstate iabeied v and The quasiparticie 
excitations in the vortex bound state have a gap by the 
hybridization of different I eigenstates as discussed beiow. 


This is a consequence of that P 2 symmetry does not pro¬ 
tect zero energy modes in a iine defect. The ampiitude 
of the excitation gap is of the order of A/Ep. Therefore, 
the gap becomes iarger when kz approach kp because the 
effective Fermi energy = (h^/2m)(kp—kl) decreases. 

Here, iet us consider the hybridization of different I 
eigenstates for non-axisymmetric v vortex. For axisym- 
metric v vortex, since anguiar momentum I is a weii- 
defined quantum number, different I eigenstates do not 
hybridize each other. The reduction of axiai symmetry 
to n-foid rotationai symmetry, however, gives hybridiza¬ 
tion between the states with anguiar momentum I and 
I -|- nm, where m G Z. Then, I and I + 2m eigenstates 
for non-axisymmetric v vortex, which has two-foid rota¬ 
tionai symmetry, are hybridized. For the fcz-dispersion of 
the originai v vortex bound state, branches of I and —I 
eigenstates are crossed on the zero energy, where I > 0 
eigenstates are shown in Fig. He) and V = —I eigenstates 
have eigenvaiues E' = —E. Since the I and —I eigenstates 
are hybridized for non-axisymmetric v vortex, the quasi¬ 
particie excitations have a gap. 

VI. w VORTEX 

The w vortex with the magnetic 7 r-rotation symmetry 
P 3 is described by the coefficients Cab(p) in Eq. (fTBl) as 
weii as V vortex. However, the coefficients C+o, Co+, 
Co-, and C-o are pure imaginary numbers. Thus, the 
reiative phase between the compensate A-phase and ,5- 
phase components and the buik B-phase component is 
7 r/ 2 , which is the difference of the OP from v vortex. 

We caicuiate the vortex bound state for w vortex by 
the BdG theory with the v vortex OP shown in Fig. 
after changing the phases of C+o, Co+, Cq-, and C_o by 
7 r/ 2 . The obtained eigenvaiues are shown in Figs. 
and|3Kd). In Fig.|3](c), the eigenvalues for fcz = 0 are clas¬ 
sified into quantized orbital angular momentum I along 
the vortex line. The vortex bound state has degenerate 
two exact zero energy modes at i = 0 indicated by an 
arrow. The fcz-dispersion of the vortex bound state for 
various Vs is shown in Fig.0](d). The quasiparticles with 
I = 0 are degenerate in small kz and linearly cross the 
zero energy at kz = 0. The I = 0 branches cross the 
zero energy twice at = 0 and \kz\ kp schematically 
shown in Fig. |D( 6 )j which is consistent with the value of 
second Ghern number discussed below. 

The wave functions of the degenerate zero energy 
modes indicated by the arrow in Fig. He) are shown in 
Figs. Ha) andHb)- The quasiparticie wave functions 

for each mode satisfy Re^uJ^ = Re^Up^ and Im^Up^ = 

—Im^UQ^, namely, uj = and Mq = . There¬ 

fore, the degenerate zero energy modes are the Majorana 
zero modes. 

The existence of the Majorana zero modes can be un¬ 
derstood by the topological arguments in Appendix |B] 
The additional magnetic 7 r-rotation symmetry P 3 in the 



















kz/kp 


FIG. 4. (Color online) Numerical results for w vortex, (a) 
and (b) Degenerate wave functions Uq{p) and Vo(p) with 
= 0. (c) Eigenvalues for quasiparticles with kz = 0. 
The degenerate eigenstates indicated by an arrow correspond 
to the two wave functions in (a) and (b). (d) fc 2 -dispersion of 
the vortex bound state whose color denotes angular momen¬ 
tum I of quasiparticles, (e) Schematic fc^-dispersion for I = 0 
quasiparticles, which is precisely degenerate near kz = 0. 

Altland-Zrinbauer symmetry class D gives a topological 
classification Z 0 Z for the vortex bound state in three- 
dimensional space (l5j |. 

One of the topological numbers Z is the second Chern 
number. For w vortex, the second Chern number is found 
to be zero as well as o vortex because the BdG Hamilto¬ 
nian for w vortex is identical to that for o vortex at infin¬ 
ity from the vortex. The vanishing second Chern number 
is consistent with the obtained fc^-dispersion of the vortex 
bound states in Fig. HKd): Some of vortex bound states 
cross the zero energy in the fc^-direction, but the cross¬ 
ing always occurs twice, so they can smoothly merge into 
the bulk state without closing bulk gap. Therefore, the 
zero energy states are topologically unstable, except for 
the I = 0 bound state. The topological stability of the 
I = 0 bound state is ensured by the other Z topological 
invariant, namely the index trF. Since the chiral sym¬ 
metry F = CV 3 is also defined in the symmetric space 
ky = kz = 0 and ^ = 0 or tt for w vortex, the one¬ 
dimensional winding numbers and are eval¬ 

uated as the same manner in Sec. El The difference of 
the winding numbers {w^ — w^^j‘1 = 2 provides the Z 
topological invariant and guarantees the existence of two 
zero energy states at = 0 , which are indeed realized 
as the two I = 0 zero modes in Fig. [H 


VII. SUMMARY 

We have calculated the bound state in o vortex, v vor¬ 
tex, and w vortex by the full quantum BdG theory with 


the self-consistent OP obtained by using the quasiclassi- 
cal theory. Moreover, we have discussed symmetry pro¬ 
tection of zero energy excitations in the vortex bound 
states with additional symmetry. Our results are sum¬ 
marized in the following and in Table L Gharacteristic 
features of the bound states in u vortex and uvw vortex 
are also listed in Table I. 

The most symmetric o vortex has Pi, P2, and P3 sym¬ 
metries in which induced components around the vortex 

core are fixed for (7++ and C _in real numbers, and 

C+o, Co+, C-o, and Gq- as zero. The vortex bound state 
for the quasiparticles with angular momentum I = 0 has 
spin degenerate Majorana zero modes at kz = 0. The 
Majorana zero modes are protected by P3 symmetry and 
characterized by a topological invariant Z which is the 
chiral index {w° — w '^)/2 = 2 . 

The V vortex has P2 symmetry in which all induced 
components around the vortex core are real numbers. In 
the vortex bound state, the quasiparticles with |Z| < 
cross the zero energy twice at finite ±/cz- However, 
the zero energy modes are topologically trivial, that 
is, they are not protected any symmetry. Then, low 
energy excitations in non-axisymmetric v vortex has a 
gap because the deformation breaking axial symmetry 
lifts the zero energy modes in spite of keeping P2 sym¬ 
metry. If we observe the excitation gap accompany¬ 
ing the vortex transition in the B-phase, that will be 
a strong evidence of that the realized vortex has only 
P2 symmetry, namely, v vortex. Although the order of 
the gap is A/Ep, that becomes large for quasiparticles 
with kz ~ kp owing to the small effective Fermi energy 
= (h^/2m)(fc| — kl). The excitation gap may have 
been observed as the difference of critical velocity be¬ 
tween VI and V2 vortices [^.1^. 

The w vortex has P3 symmetry. The difference be¬ 
tween w vortex and v vortex is only phases of induced 
components G+o, G 0 - 1 -, G_o, and Gq- concerning the OP. 
However, the vortex bound state in superfluids belonging 
to the symmetry class D with additional symmetry P 3 is 
characterized by the topological invariants Z 0 Z, the sec¬ 
ond Chern number and the chiral index. For w vortex, 
the second Chern number Ch 2 = 0, but the chiral in¬ 
dex (?ii° — w'^)/2 = 2 indicating difference between the 
numbers of the zero energy mode with opposite chirality. 
Thus, the two-fold degenerate Majorana zero modes at 
the symmetric point kz = 0 and I = 0 are protected by 
P 3 symmetry. 

The uvw vortex without any Pi, P 2 , and P 3 sym¬ 
metries has many zero energy modes, as calculated by 
Silaev [S^. However, the zero energy modes are topolog¬ 
ically unstable because the vortex bound state in super¬ 
fluids belonging to the class D without additional symme¬ 
try is characterized by the second Chern number, which 
is zero for the vortex bound state in the superfluid ^He 
B-phase. The u vortex with Pi symmetry is achieved 
by adding imaginary parts to real G++ and G_compo¬ 

nents in o vortex. The imaginary parts lift the Majorana 
zero modes at the symmetric point kz = 0 and I = 0 
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TABLE I. Classified vortices by Pi, P 2 , and P 3 symmetries which fix the element of induced components 
Cabip)- For each vortex, we show the presence or absence of zero energy modes (ZEM) and Majorana zero 
modes, as well as topological invariant classifying the vortex bound state if it is present. The candidates 
for the topological invariant are the chiral index N = (w^ — w'^)l2 = 2 and the second Chern number 
Ch2 = 0. 


vortex 

symmetry 

C++, C— 

C+o, Co+, 
C-o, Co- 

ZEM 

Majorana 

top. inv. 

0 vortex 

Pi, P 2 , P 3 

Real 

- 

/ 

/ 

N 

axisym. v vortex 

P 2 

Real 

Real 

/ 

- 

- 

non-axisym. v vortex 

P 2 

Real 

Real 

- 

- 

- 

w vortex 

P 3 

Real 

Imaginary 

/ 

/ 

iV©Ch 2 

uvw vortex 

- 

Complex 

Complex 

/ 

- 

Ch 2 

u vortex 

Pi 

Complex 

- 

/ 

- 

- 


by breaking P 3 symmetry. Other accidental zero modes 
will be left but they are not characterized by the second 
Chern number which generally vanishes under Pi sym¬ 
metry. 

In conclusion, Majorana zero modes protected by P 3 
symmetry are bound in o vortex and w vortex which are 
not stable in the bulk B-phase. The v vortex is more sta¬ 
ble than the o vortex owing to the condensation energy 
by A-phase and /3-phase components which compensates 
the vortex core. Then, when we confine the B-phase in 
a thin slab with hight along the vortex line to suppress 
the /3-phase component and simultaneously apply a mag¬ 
netic field along the vortex line to suppress the A-phase 
component, o vortex with Majorana zero modes will be 
achieved. Note that the confinement and the magnetic 
held do not break any Pi, P 2 , and P 3 symmetries. 
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sectors and each sector has a Majorana zero mode. The 
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Appendix A: Linear fc^-dispersion of o vortex bound 

state 


We show the linear fc^-dispersion of the vortex bound 
state in o vortex on the basis of perturbation theory [ 3 ^ . 


where uj(p) = A/'-t-Jo(^Fp) exp(—p/^ —f 7 r/ 4 ) and Uq{p) = 
■N'iJi{kpp) exp(—p/^ -|- f 7 r/ 4 ) with normalization factors 
Ma and the n-th order Bessel function Jn [11,14^. The 
wave function of the Majorana zero mode in each spin 
sector differs in the phase winding and the order of the 
Bessel function since the Hamiltonian for the f-spin (),- 
spin) sector corresponds to the Hamiltonian for the chiral 
p-wave superfluid with the chirality antiparallel (parallel) 
to the vorticity. 

Here, we consider the /c^-dispersion with small kz <?/ kp 
by perturbation theory. The perturbation Hamiltonian 
is described by 


■ffp(p) = -^Ao(p) 


/ 0 
0 
0 


0 0 P‘>‘\ 

0 0 

e-*'^ 0 0 ’ 

0 0 0 / 


(A4) 
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where projections of Hp onto Mg and Ug give 


/ 

4(p) 

V 

4(p) 


Hpip)uUp) 

^p(p)4(p) 
k. 


4(p) 

4(p) 


Hp{p)ul{p) 


Hp{p)ul{p) 


= 2A/'tA/j,^Ag(p)Jg(fcFp)^i(fcF/o)e (A5) 


Then, the fc^-dispersion, 

E± = ± 47 rA/'tA/'^^ J dpAo{p)Jo{kFp)Ji{kFp)t 


-2p/« 


(A 6 ) 


is linear to kz ■ The 0-th order perturbed wave functions 
are 

Uq (p) = -T [4 (p) ± 4(p)] . ( A7) 

These wave functions u^{p) also indicate self-conjugate 
Majorana quasiparticles. 


Appendix B: Topological classification of vortex 
bonnd states by Clifford algebras 


HBdGik, 4>) = VdHBAG[-kx,ky, k ^,\ where Va = 
i&zTzIC. Among these discrete symmetries, we respect 
P 2 and/or P 3 in the following, since Pi does not provide 
symmetry protected defect zero modes in general. Below 
we identify all topological invariants relevant to existing 
vortex zero modes in the B-phase. 

First note that the BdG Hamiltonian far way from 
the vortex core can be written in terms of the gamma 
matrices 7 ^ = {-azT^, -Ty, -Tz) as 

HBdG(fe, </) = - ^ 74 , (B3) 

with dy, = {ABkx/kF,ABky/kF,ABkz/kF,h'^kp/2m). 
The gamma matrices 7 ^ obey the Clifford algebra, 
{ 7 p: 7 i/} = ‘2-5y,v Whereas the BdG Hamiltonian near 
the core is different from Eq. (IB31) . it is smoothly inter¬ 
polated from this, with keeping symmetry of the vor¬ 
tex. Therefore, the topological structure of a vortex can 
be evaluated from the asymptotic Hamiltonian Ea. (IB3l) . 
subject to a set of symmetries of the vortex. 

To clarify the possible topological structure, we fur¬ 
thermore deform Eq. (IB3I) into the form of a Dirac Hamil¬ 
tonian. For this purpose, we may add the following topo¬ 
logically trivial system with the same discrete symme¬ 
tries, 


In this appendix, we discuss topological classification 
of vortex bound states in the B-phase. Consider a vortex 
along the z-axis, and a circle surrounding the vortex, 
which is parametrized by the angle evaluated from the 
x-axis. The semi-classical BdG Hamiltonian on the circle 
is given by 


H-BdGik, (j)) 


( h{k) A{k,4i) \ 


(Bl) 


where h{k) = {It? /2m){k^ — fc|)i is the Hamiltonian in 
the normal state, and A is OP which approaches 


A(fc,</) 


Ab 

fcp 


-kx + iky kz 

kz kx + iky 


P*, (B2) 


far away from the vortex core. The BdG Hamiltonian has 
particle-hole symmetry defined by CHBdG{k,4>)C~^ = 
—HBdG{—k,(j)) with C = TxJC and = 1, where K. is 
the complex conjugation operator. 

The axisymmetric vortices in the ^He B-phase may 
have three types of discrete symmetries: Pi, P 2 , 
and P 3 [ 1 ^ . The inversion symmetry Pi implies 
PBdG(fc,</) = PBdG(-fe,/'+ tt). P 2 is magnetic re¬ 
flection symmetry that obtained as a combination of 
time-reversal and mirror reflection with respect to a 
plane including the vortex line: If we take the xz- 
plane as the reflection plane, P 2 reads HBdG{k,(j)) = 
P 2 HBdG{-kx,ky,-kz,-(l))p 2 ^^ where P 2 = dyTzf = 
iTzK. with the time-reversal operator P = idyJC. P 3 is 
magnetic 7r-rotation symmetry around an axis perpen¬ 
dicular to the vortex line, say, the i-axis. P 3 means 




g-d0/2)74^^...y^gi(<^'/2)74 




2 m 


74- 


(B4) 


In comparison with the original BdG Hamiltonian, 
i7gdG(^’ (/) has a negative chemical potential —h?kp/2m, 
and thus it is deformable to a topologically trivial in¬ 
sulator without gap closing, by taking the limit Ab —>■ 
0. This means that PggQ(fc,(/) is topologically triv¬ 
ial. Adding the topological trivial band to the orig¬ 
inal one, we have the extended BdG Hamiltonian 
HeBdG{k, 4>) given by the direct product of PBdG(fc;</) 
and H^^Qik,(j)), 


HeBdG{k, 4>) 


HBdG{k,(j)) ^ 0 


, (B5) 


which is again written in terms of the gamma matrices 
as 


PeBdG(fc,/>) = 


H^k^ 

2 m 


74, (B 6 ) 


where the new gamma matrices F^ are given by F^ = 
ItiPz with the Pauli matrix pi in the band space. The 
particle-hole and other discrete symmetries imply that 

[c, Fi] = [c, r 2 ] = [c, F 3 ] = {c, r 4 } = o, (B7) 
{P 2 ,ri} = [P 2 ,r 2 ] = {P2,r3} = [P2,r4] = o, (bs) 
{P 3 ,ri} = [P3,r2] = [P3,r3] = [P3,r4] = o, (B9) 

with 


[C,P2] = [C,P3] = [P2,P3] = 0, 


(BIO) 
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(74 obeys the same (anti-)commutation relations as r 4 .) 
The extended BdG Hamiltonian is stable-equivalent to 
the original one in the sense of the K-theory. 

For the extended BdG Hamiltonian, we can introduce 
Fs = that satisfies 

Tl = i, {r5,r4 = o, (Bii) 

{C, Fs} = {r 2 , F 5 } = {Vs, Fs} = 0 . (B 12 ) 

Using F 5 , we perform the deformation of the Hamiltonian 
that preserves all the discrete symmetries (as well as the 
particle-hole symmetry) except for Pi: 

'' h'^k'^ 

H^{k,4,) = C/a(<(>)d^F^C/-^(</>) - ^ 74 , (B13) 

where Ua{ 4 >) is given by 

Ua{4>) = e®(';^/2)[cosO!r4+isinar4r5]_ (B14) 


Through this equation, the extended BdG Hamiltonian 
at a = 0 is smoothly deformed into the following Dirac 
Hamiltonian at a = 7 r/ 2 , 

iJD(fc, 4>) = kxVi + kyT2 + kj'^ + cos(()r4 -I- sinc/iFs, 

(B15) 


where we have placed Ab/A:f = 1 and c ?4 = 1 and omit¬ 
ted the regularization term —{h^k^/ 2 m)"fi for simplicity. 
The Dirac Hamiltonian has the same topological proper¬ 
ties as the original BdG Hamiltonian. 

To elucidate the topological structure of the Dirac 
Hamiltonian, we consider a family of Dirac Hamiltoni¬ 
ans which have the same form and the same symmetries 
as Eq. (IB15|) . Even for these Dirac Hamiltonians, F^, 
Fs, C, 7 ^ 2 , and V 3 should satisfy the same commutation 
or anti-commutation relations as Eqs. (IB7I) - (IB10|) . but 
their matrix representation are not specified anymore. 
For the family of Dirac Hamiltonians, F 4 transforms like 
a mass term under the symmetries, so we can clarify the 
topological structure by counting topologically distinct 
r 4 matrices consistent with the symmetries |14.1^. 

To count topologically distinct r 4 matrices, we use the 
fact that the symmetry operators and gamma matrices 
form a real Glifford algebra Clp^q that hasp-fg generators 
{ei,..., Cp; Cp+i,..., ep+q} satisfying 

{ci, Cj} = 0 , {i ^ j), 


f - 1 , 1 < i < p 

( -1-1, p-|-l<i<p-|-g 


(B16) 


For instance, consider uvw vortex that only has the 
particle-hole symmetry. The particle-hole symmetry C 
and the gamma matrices other than r 4 form as 


{TiJ,T 2 J,T 3 J-,C.CJ,T 3 }, (B17) 


with 

= c^ = i, (cjf = i, rl = i. (BI 8 ) 


Here we have introduced representing the imaginary 
unit “i” so that we can treat complex structure originated 
from the anti-unitary operator C. If we take into account 
r4 as well, the Clifford algebra is extended to 

{ri, 7 ,r 2 , 7 ,r 3 , 7 ;C,c, 7 ,r 5 ,r 4 }. (bi 9 ) 

Therefore, having a r4 matrix consistent with the 
particle-hole symmetry provides an extension of the Clif¬ 
ford algebra from C/3,3 to C/3,4, and vice versa. A set of 
the latter extensions defines the classifying space TZq, so 
topologically distinct r4 matrices can be counted as the 
number of the disconnected parts of the classifying space, 
i.e. 7 ro( 7 ^o) = Correspondingly, we can introduce 
the second Chern number defined in Eq. ^ HI, [ 43 . 
For a vortex in ^He-B phase, however, the second Chern 
number becomes zero since its asymptotic Hamiltonian 
Eq. (IB 3 I) has an additional inversion symmetry. (See also 
discussions in Sec. El) 

For o vortex, the corresponding Dirac Hamiltonian 
should be subject to additional P 2 and P3 symmetries. 
In the presence of P 2 and P3, symmetry operators and 
the gamma matrices other than r4 form C/4,4 as 

{TiJ, T 2 J, T 3 J,CT2V3V2-,C,CJ, T3,JT3V2V3}, 

(B 20 ) 

then a r4 matrix extends this into C/4,5 as 

{TiJ, T 2 J, T 3 J, ct2T3V2\C,cj, T3,jt3V2V3, r 4 }. 

(B 2 I) 

As the classifying space for the extension C/4,4 —>■ C/4,5 
is i?o, the BdG Hamiltonian with o vortex is classified as 
7>'o(-/?o) = The correspondin g to pological number is 
(u>° — w^) 12 defined in Eq. (1^ jl5j|. 

For V vortex, the Dirac Hamiltonian has additional P 2 
symmetry. Possible r4 matrices can be identified as the 
extension of Clifford algebra C/4,3 —>■ C/4,4: 

{TiJ, T 2 J, T 3 J, CT2T3V2\C, CJ, r 5 } 

^ {TiJ,T 2 J,T 3 J.CT 2 T 3 V 2 -,C,CJ,T 3 ,Ti}. (B 22 ) 

The classifying space for the extension C/4,3 C/4,4 is 

i?_i ~ i?7. Since 7ro(i?7) = 0 , u vortex does not support 
a zero mode protected by P2 symmetry. 

Finally, consider w vortex. In this case, the additional 
symmetry is P3. Topologically distinct r4 matrices can 
be identified by using the extension 

{TiJ, T 2 J, T,J- C, CJ, r5} O {V' 3 } 

^ {TiJ, T 2 J, T,J-C,CJ, Fs, r4} 0 (B 23 ) 

with P3 = CJT2r3r5P3. This gives C/3,3 0 C/0,1 —>■ 
C/3,4 0 C/0,1, which classifying space is given by i?o x 
i?o- Therefore the topological classification of the BdG 
Hamiltonian with w vortex is ttq(Rq x Rq) = Z 0 Z. 
The corresponding topological invariants are the second 
Ghern number in Eq. ( 1271 ) and the one-dimensional wind¬ 
ing number defined in Eq. (ESI). 
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Appendix C: Phase shift across a vortex core 


In this appendix, we make a detailed explanation for 
the excitation spectrum shown in Fig. EJe) in terms of 
the OP structure. When a vortex core is compensated by 
superfluid components, the vortex bound state for quasi¬ 
particles across the vortex core can be regarded as the 
Andreev bound state on the junction with the compen¬ 
sated superfluid between the bulk B-phases. If the OPs 
of them are simultaneously diagonalizable on a quasipar- 

I 


tide path, we can estimate the bound state energy from 
the phase shift on the vortex core. 

The explicit OPs in the bulk B-phase and for the in¬ 
duced superfluid in the vicinity of an axisymmetric vortex 
are described by 


AB(fc,p) = AB(p)e*‘^ 


svaOke 
cos 9 k 


cos 6k 
sin 


(Cl) 


and 


A C'-so(p)cos6»fc-tC'++(p)sin6»fee*('^‘^ 

core! \^Co+(p)sin6»fce*'^"-C'o-(p)sin6»fce*(2<^-’^'= 


where (j)k and 6k denote the azimuthal and polar an¬ 
gles of the quasiparticle momenta, respectively, and (j) 
is an azimuthal angle of the real space coordinate. Since 
slightly induced components C++ and C_do not in¬ 

fluence phase shifts of the OP, we disregard them in the 
following discussion. 

Let us focus our attention on the quasiparticle path 
across the vortex core shown in Fig.[5](b). The azimuthal 
angles of the quasiparticle momentum and the real space 
are 4> = (j)k or (f> = (j)k + tt, where the quasiparticle has 
the angular momentum I = 0. For the quasiparticles with 
1 = 0, the OP of the vortex core state in Eq. (E 2 l) is 
simplified as 

A^,,{k,p) = AA(0fc,p)d^'“<T, + Ap{6k,p)^^ 

-bA_o(0fc,p)d'^''e-*^'“"O (C3) 

where Ajk{6k,p) = (Co+ - Co-)sin 6 »fc, Ap{6k,p) = 
(C+o - C_o)cos 6 >fe, and A_o( 6 'fc,p) = C_ocos 6 »fc. Fig¬ 
ure [SKa) shows the 6k- and p-dependence of the OP com¬ 
ponents Ajk{6k,p), Ap{6k,p), and A_o( 6 >fc,p), which are 
obtained from our numerical results in Fig. EJa). 

As shown in Fig. [SKa), the A-phase component A a 
dominates at the vortex core for 0fc = 7 r/ 2 . The quasipar¬ 
ticles with momentum perpendicular to the vortex line 
can be regarded as quasiparticles across the junction with 
the A-phase between the bulk B-phase. The OPs of the 
bulk B-phase in Eq. (EU and the A-phase described by 
the A+ term in Eq. (IC3I) are simultaneously diagonaliz¬ 
able by using the unitary matrix 

t/ = exp [—i( 7 r/ 4 )(fc • <t)] . (04) 

Then, they are transformed as 

UA^{(l)k,6k = 7r/2,p, (/) = -I- 

= ±AB(p)C^'“e-*^^"^ ly (05) 


Co+(p)sin0fcC‘^''-Co-(p)sin0feC(2<^-'^'=) \ . . 

) C_o(p)cos6»fee2fo-C__(p)sin6»fee*(3<^^-‘^)y ’ ^ 



20 0.0 


FIG. 5. (Color online) (a) OP components Aa (dash-dotted 
line), Ap (solid line), and A_o (dashed line) in Eq. (IC3II as a 
function of the quasiparticle momentum 9k and distance from 
a vortex core p. (b) Thick arrow denotes a quasiparticle path 
defined by 9k and (j>k- 


where -|- (—) sign is taken for (j) = (j)k {(t> = (t>k + tt), and 

C/AA( 0 fe = ^/ 2 ,p)e*-^'=a,C/^ 

= -[Co+(p) - ( 0 * ) ’ 

Eor V vortex with real Cqa ^^^d Cq-, comparison of 
Eqs. (|05ll and (1061) shows that the quasiparticles feel 
the phase shift p = ± 7 r /2 for each spin sector. The 
eigenenergy of the Andreev bound states on a junction 
with phase shift p is E oc ±| cos((p/2)| Therefore, 

quasiparticles with momentum dk = 7r/2 and angular 
momentum I = 0 have an energy gap in the vortex bound 
state. This gap corresponds to the energy gap at fcz = 0 
and / = 0 in Figs.EKd) andEIe). 
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On the other hand, Fig. [SJa) shows that the A_o term 
dominates around the vortex core for the quasiparticle 
with momentum 0^ = 0. Note that A_o and A^ induced 
on the vortex core have the same phase factor. We can 
simultaneously diagonalize the OPs of the bulk B-phase 
in Eq. den and the A_o term in Eq. (jC3p by using the 
unitary matrix 


U = exp 






(C7) 


where = —aisini^fc + ycos(j)k- After the unitary 
transformation at dk = 0, the OP of the A_o term re¬ 
mains unchanged and that of the bulk B-phase trans¬ 
forms into the same form in Eq. (IC5I) . Therefore, quasi¬ 
particles with momentum k almost parallel to the vortex 
line feel the phase shift tp = tt. This 7r-phase shift results 
in the I = 0 excitation spectrum approaching the zero 
energy at = i/cp as shown in Eig. ISIe). 

In this sense, the A_o term plays a crucial role in the 


excitation spectrum. If we set only C+o a-nd Co+ compo- 
nents to be non-zero near the vortex core as Ref. [30j, 
we can diagonalize the OP of the bulk B-phase and 
compensated components in the vortex core simulta¬ 
neously by using the unitary matrix in Eq. (IC7)) with 
9k = arccos[±-\/2Co+/(2C'o+ — C+o)]: where (70+ and 
C+o have opposite signs. Since quasiparticles with the 
momenta = ife+lcos^fcl feel the 7r-phase shift, the 
I = 0 excitation spectrum should cross the zero energy in 
0 < Ifc^l < fcp- However, we can obtain I = 0 excitation 
spectrum crossing the zero energy at kz = ±fcp because 
we consider not only C+o and Co+ but also the other 
induced components. 

The phase shift also explains the existence of the de¬ 
generate Majorana zero modes in w vortex. For w vor¬ 
tex, since Co+ and Cq- in Eq. (ICel) are pure imaginaries, 
quasiparticles with kz = 0 and I = 0 feel the 7r-phase 
shift, which is confirmed by comparison of Eqs. (IC5I) and 
(IC6D . The zero energy bound states at the particle-hole 
symmetric point owing to the 7r-phase shift are exactly 
degenerate Majorana zero modes. 
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